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We study the supersymmetry reach of the Tevatron in channels containing both isolated lep- 
tons and identified tau jets. In the most challenging case, where the branching ratios of gauginos 
to taus dominate, we find that searches for two leptons, a tau jet and a large amount of missing 
transverse energy have a much better reach than the classic trilepton signature. With total inte- 
grated luminosity of L > 4 fb _1 , the Tevatron will start extending the expected LEP-II reach for 
supersymmetry. 
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Q\ ' Searches for supersymmetry (SUSY) in Run I of the Tevatron have been done exclusively in channels involving 
ON ■ some combination of leptons, jets, photons and missing transverse energy (ffr) 0- At the same time, several Run I 
analyses have identified hadronic tau jets, e.g. in W-production and top decays j|. Hadronic taus have also been 
used to place limits on a charged Higgs Q and leptoquarks ||. As tau identification is expected to improve further 
Qh' in Run II, this raises the question whether SUSY searches in channels involving tau jets are feasible. 
^ , SUSY signatures with tau leptons are very well motivated, since they arise in a variety of models of low-energy 
supersymmetry, e.g. gravity- mediated |3p|| or the minimal gauge-mediated models [||-@. In this paper we shall 
study all possible experimental signatures with three identified objects (leptons or tau jets) plus J^t-, and compare 
their reach to the clean trilepton channel jll],[l2| , which is the classic SUSY signature at the Tevatron. It arises in the 
, decays of gaugino-like chargino-neutralino pairs xfx2- The reach is somewhat limited by their rather small leptonic 
' branching fractions. In the limit of either heavy or equal in mass squarks and sleptons, the leptonic branching ratios of 
Xi and %2 are W'-like and Z-like, respectively. However, both gravity-mediated and gauge-mediated models of SUSY 
breaking allow the sleptons to be much lighter than the squarks, thus enhancing the leptonic branching fractions of 
, the gauginos. 

There are at least three generic reasons as to why one may expect light sleptons in the spectrum. First, the slepton 
masses at the high-energy scale may be rather small to begin with. This is typical for gauge-mediated models, since the 
Q\ ' sleptons are colorless and do not receive large soft mass contributions proportional to the strong coupling constant 
a s . The minimal gravity- mediated (mSUGRA) models, on the other hand, predict light sleptons if the universal 
O |' scalar mass M is much smaller than the universal gaugino mass M 1 / 2 - Second, the renormalization group equations 
for the scalar soft masses contain terms proportional to Yukawa couplings, which tend to reduce the corresponding 
mass during the evolution down to low-energy scales. This effect is significant for third generation scalars, and for 
large values of tan/3 (the ratio of the Higgs vacuum expectation values v 2 and v±) splits the staus from the first 
two generation sleptons. And finally, the mixing in the charged slepton mass matrix further reduces the mass of the 
lightest eigenstate. The slepton mixing is enhanced at large tan/3, since it is proportional to [imi tan /3/m~, where mi 
(mj) is the lepton (slepton) mass and \i is the supersymmctric Higgs mass parameter. Notice that this effect again 
£3 , only applies to the staus, since m T ;§> m M e . 

Due to these three effects, it may very well be that among all scalars, only the lightest sleptons from each generation 
(or just the lightest stau fi) are lighter than xt an d xH- Indeed, in both gravity-mediated and gauge- mediated models 
one readily finds regions of parameter space where either m^o < rrif 1 ~ mp, R < rn^+ ~ to^o (typically at small tan /3) 
or to^o < rrif 1 < rn^+ ~ rn^o < m^ R (at large tan /?) . Depending on the particular model, and the values of the 
parameters, the gaugino pair decay chain may then end up overwhelmingly in any one of the four final states: III, Ut, 
Itt or TTT. 

In order to make a final decision as to which experimental signatures are most promising, we have to factor in the 
tau branching ratios to leptons and jets. About two-thirds of the subsequent tau decays are hadronic, so it appears 
advantageous to consider signatures with tau jets in the final state as alternatives to the clean trilepton signal. (From 
now on, we shall use the following terminology: a "lepton" (I) is either a muon or an electron; a tau is a tau-lepton, 
which can later decay either leptonically, or to a hadronic tau jet, which we denote by t/j.) The presence of taus in 
the underlying SUSY signal always leads to an enhancement of the signatures with tau jets in comparison to the clean 
trileptons. This disparity is most striking for the case of rrr decays, where BR{ttt — > Ut^) / BR(ttt — > III) ~ 5.5. 
An additional advantage of the tau jet channels is that the leptons from tau decays are much softer than the tau jets 
and as a result will have a relatively low reconstruction efficiency. 

On the other hand, the tau jet channels suffer from larger backgrounds than the clean trileptons. The physical 
background (from real tau jets in the event) is actually smaller, but a significant part of the background is due to 
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events containing narrow isolated QCD jets with the correct track multiplicity, which can be misidentified as taus. 
The jetty signatures are also hurt by the lower detector efficiency for tau jets than for leptons. The main goal of our 
study, therefore, will be to see what is the net effect of all these factors, on a channel by channel basis. 

For our analysis we choose to examine one of the most challenging scenarios for SUSY discovery at the Tevatron. 
We shall assume the typical large tan/3 mass hierarchy to^o < < < mji R - One then finds that BR(x_iX2 ~* 
ttt + X) ~ 100% below xf — > W^Xi an( A X2 ~~ * %Xi thresholds. In order to shy away from specific model 
dependence, we shall conservatively ignore all SUSY production channels other than X\Xj2 P a i r production. The pt 
spectrum of the taus resulting from the chargino and neutralino decays depends on the mass differences m -+ — 
and mfj — m^o. The larger they are, the harder the spectrum, and the better the detector efficiency. However, 

as the mass difference gets large, the xt anc ^ X2 masses themselves become large too, so the production cross- 
section is severely suppressed. Therefore, at the Tevatron we can only explore regions with favorable mass ratios 
and at the same time small enough gaugino masses. This suggests a choice of SUSY mass ratios: for definiteness 
we fix 2m^o ~ (4/3) ~ m~+(< mfi B ) throughout the analysis, and vary the chargino mass. The rest of the 
superpartners have fixed large masses corresponding to the mSUGRA point M = 180 GeV, M 1 / 2 — 180 GeV, A = 
GeV, tan j3 = 44 and fi > 0, but we are not constrained to mSUGRA models only. Our analysis will apply equally to 
gauge-mediated models with a long-lived neutralino NLSP, as long as the relevant gaugino and slepton mass relations 
are similar. Note that our choice of heavy first two generation sleptons is very conservative. A more judicious choice 
of their masses, namely mf lR < would lead to a larger fraction of trilepton events, and as a result, a higher 

reach. Furthermore, the gauginos would then decay via two-body modes to first generation sleptons, and the resulting 
lepton spectrum would be much harder, leading to a higher lepton efficiency. Notice also that the X1X2 production 
cross-section is sensitive to the squark masses, but since this is the only production process we are considering, our 
results can be trivially rescaled to account for a different choice of squark masses, or to include other production 
processes as well. 

Since the experimental signatures in our analysis contain only soft leptons and tau jets, an important issue is 
whether one can develop efficient combinations of Level 1 and Level 2 triggers to accumulate these data sets without 
squandering all of the available bandwidth. We will not attempt to address this issue here; instead we will assume 
100% trigger efficiency for those signal events which pass all of our analysis and acceptance cuts. We have nevertheless 
studied the following set of triggers: 1) # T > 40 GeV; 2) p T (l) > 20 GeV and 3) p T (l) > 10 GeV, p T (jet) > 15 GeV 
and-^r > 15 GeV; with pseudorapidity cuts |r?(e)| < 2.0, \f](fj,)\ < 1.5 and \n(jet)\ < 4.0. We found that they are 
efficient in picking out about 90 % of the signal events in the channels with at least one lepton (see below) . Dedicated 
low pt tau triggers for Run II, which may be suitable for the new tau jet channels, are now being considered by both 
CDF H| and DO |l|. 

We used PYTHIA v6.115 and TAUOLA v2.5 for event generation. We used the SHW v2.2 package @, which 
simulates an average of the CDF and DO Run II detector performance. In SHW tau objects are defined as jets with 
\rj\ < 1.5, net charge ±1, one or three tracks in a 10° cone with no additional tracks in a 30° cone, Et > 5 GeV, 
Pt > 5 GeV, plus an electron rejection cut. SHW electrons are required to have \rj\ < 1.5, Et > 5 GeV, hadronic to 
electromagnetic energy deposit ratio Rh/ e < 0.125, and satisfy standard isolation cuts. Muon objects are required to 
have \r]\ < 1.5, Et > 3 GeV and are reconstructed using Run I efficiencies. We use standard isolation cuts for muons 
as well. Jets are required to have |?7| < 4, Et > 15 GeV. In addition we have added jet energy correction for muons 
and the rather loose id requirement Rh/ e > 0.1. We have also modified the TAUOLA program in order to correctly 
account for the chirality of tau leptons coming from SUSY decays. 

The reconstruction algorithms in SHW already include some basic cuts, so we can define a reconstruction efficiency 
e rec for the various types of objects: electrons, muons, tau jets etc. We find that as we vary the chargino mass from 
100 to 140 GeV the lepton and tau jet reconstruction efficiencies for the signal range from 42 to 49 %, and from 29 
to 36%, correspondingly. The lepton efficiency may seem surprisingly low, but this is because a lot of the leptons are 
very soft and fail the Et cut. The tau efficiency is in good agreement with the results from Ref. |l6| and jL7), once 
we account for the different environment, as well as cuts used in those analyses. 

The most important background issue in the new tau channels is the fake tau rate. Several experimental analyses 
try to estimate it using Run I data. Here we simulate the corresponding backgrounds to our signal and use SHW 
to obtain the fake rate, thus avoiding trigger bias |l6|]. We find that the tau fake rate in W production is 1.5%, 
independent of the tau pt, which is in agreement with the findings of Refs. JToj [l8t] - 

In the following we list our cuts for each channel. 

In order to maximize the reach in the IUE/t channel, we apply the soft lepton pt cuts advertised in Refs. |ll}| . We 
require a central lepton with pt > 11 GeV and < 1.0, and in addition two more leptons with pt > 7 GeV and 
Pt > 5 GeV. Leptons have to be isolated: 1(1) < 2 GeV, where / is the total transverse energy contained in a cone of 
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size SR = \J Aip 2 + Ar] 2 = 0.4 around the lepton. We impose a dilepton invariant mass cut for same flavor, opposite 
sign leptons: \mu — Mz\ > 10 GeV and \mu \ > 11. Finally, we impose an optional veto on additional jets and require 
1$t to be either more than 20 GeV, or 25 GeV. This gives us a total of four combinations of the T$t cut and the jet 
veto (JV) (A:#r > 20 GeV, no JV; B:^ T > 25 GeV, no JV; C:^ T > 20 GeV, with JV; D:^ T > 25 GeV, with JV), 
which we apply for all tau jet signatures later as well. 

For our Uth $t analysis we impose cuts similar to the stop search analysis in the jtfr channel pj[ : two isolated 
(/(/) < 2 GeV) leptons with pt > 8 GeV and px > 5 GeV, and one identified tau jet with pt{t~}i) > 15 GeV. Again, we 
impose the above invariant mass cuts for any same flavor, opposite sign dilepton pair. This channel was also studied 
in Ref. |?J with somewhat harder cuts on the leptons. A separate, very interesting signature (Z + Z + r^^r) arises if the 
two leptons have the same sign, since the background is greatly suppressed. In fact, we expect this background to be 
significantly smaller than the trilepton background! Roughly one third of the signal events in the general IWh sample 
are expected to have like-sign leptons. 

For our IthTk $t analysis we use some basic identification cuts: two tau jets with pt > 15 GeV and px > 10 GeV 
and one isolated lepton with pt > 7 GeV. 

Finally, for the t^t^t^ ]£t signature we only require three tau jets with px > 15, 10 and 8 GeV, respectively. 

One can get a good idea of the relative importance of the different channels by looking at the corresponding signal 
samples after the analysis cuts have been applied. In Fig. [l] we show the signal cross-sections times the corresponding 
branching ratios times the total efficiency e to t = ^rec^cuts, which accounts for both the detector acceptance e rec and 
the efficiency of the cuts e C uts (for each signal point we generated 10 5 events). We see that the lines are roughly 
ordered according to the branching ratios of three taus into the corresponding final state signatures. This can be 
understood as follows. The acceptance (which includes the basic ID cuts in SHW) is higher for leptons than for r jets. 
Therefore, replacing a lepton with a tau jet in the experimental signature costs us a factor of ~ 1.5 in acceptance, due 
to the poorer reconstruction of tau jets, compared to leptons. Later, however, the cuts tend to reduce the leptonic 
signal more than the tau jet signal, mostly because the leptons are softer than the tau jets. It turns out that these 
two effects mostly cancel each other, and the total efficiency etot is roughly the same for all channels. Therefore the 
relative importance of each channel will only depend on the tau branching ratios and the backgrounds. For example, 
in going from III to Utu, one wins a factor of 5.5 from the branching ratio. Therefore the background to IWhl^T must 
be at least 5.5 2 ~ 30 times larger in order for the clean trilepton channel to be still preferred. 

We next turn to the discussion of the backgrounds involved. We have simulated the following physics background 
processes: ZZ, WZ, WW, ti, Z+jets, and W+jets, generating 4 x (10 6 ) and 2 x (10 7 ) events, respectively. We 
list the results in Table |, where we show the total background cross-section oqq for each case A-D, as well as the 
contributions from the individual processes, for case A. All errors are purely statistical. Our simulated background 
in the trilepton channel is higher than previously found in Refs. [0,0, which employed ISAJET for event generation. 
Current versions of ISAJET simulate the WZ and ZZ processes in the limit of zero Z width. We find that most of 
the W Z and ZZ background events from PYTHIA contain an off-shell Z instead and pass the dilepton invariant mass 
cuts. Unfortunately, neither ISAJET, nor PYTHIA contain the Wj interference contribution to WZ, so our result 
still somewhat underestimates the WZ trilepton background. As we move to the channels with tau jets, the number 
of events with real tau jets decreases, mostly because of the smaller branching ratios of W and Z to taus. However, 
the contribution from events with fake tau jets increases, and for the 3t/j channel events with fakes are the dominant 
part of the background. Our results for the Ut^ and It^t^ channels differ from those of Ref. (To), which assumed a 
p<r-dependent fake tau rate of only 0.1-0.5%. 

We find that although the jet veto successfully removes the ti background to the first three channels, it also reduces 
the signal (see Fig. [l]), and hence does not improve signal-to- noise. A higher fyr cut also never seemed to help. Indeed, 
the major backgrounds contain leptonic t and W decays and tend to have a lot of missing energy. Notice that we did 
not account for fake leptons in our Z + jets backgrounds to 3/ and l + l + Th, since SHW does not provide a realistic 
simulation of those. The best way to estimate this background component will be from Run II data (an analysis 
based on Run I data [pOfl reveals that the 3^ contribution from fakes is comparable to our result in Table [|, which only 
includes real isolated leptons from heavy flavor jets). We have a bit underestimated the total background to the 3r^ 
channel by considering only processes with at least one real tau in the event. We expect sizable contributions from 
pure QCD multijet events, or Wj — ► jjj, where all three tau jets are fake. 

A 3cr exclusion limit requires a total integrated luminosity L(3a) = 9<jbg {"sig BR{XiX2 ~~ * ^0 £ tot} ■ Notice 
that L(3cr) depends linearly on the background obg after cuts, but quadratically on the signal branching ratios. This 
allows the jetty channels to compete very successfully with the clean trilepton signature, whose branching ratio is quite 
small. In Fig. ^ we show the Tevatron reach in the three best channels: trileptons (x), dileptons plus a tau jet (□) 
and like-sign dileptons plus a tau jet (O). We see that the two channels with tau jets have a much better sensitivity 
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compared to the usual trilepton signature. Assuming that efficient triggers can be implemented, the Tevatron reach 
will start exceeding LEP II limits as soon as Run II is completed and the two collaborations have collected a total 
of 4 ftU 1 of data. Considering the intrinsic difficulty of the SUSY scenario we are contemplating, the mass reach 
for Run III is quite impressive. One should also keep in mind that we did not attempt to optimize our cuts for the 
new channels. For example, one could use angular correlation cuts to suppress Drell-Yan, transverse W mass cut to 
suppress WZ pi] ], or (chargino) mass-dependent pr cuts for the leptons and tau jets, to squeeze out some extra reach. 
In addition, the Ut^ channel can be explored at smaller values of tan f3 as well [[z| Jl~T|j20[ ] , since the two-body chargino 
decays arc preferentially to tau sleptons. In that case, the clean trilepton channel still offers the best reach, and a 
signal can be observed already in Run II. Then, the tau channels will not only provide an important confirmation, 
but also hint towards some probable values of the SUSY model parameters. 
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TABLE I. Total background cross-section after cuts <j B g (in fb) for the various channels in cases A-D, as well as the 
breakdown of the individual contributions for case A. 



&BG (fb) 


UW/t 


llriffr 


l + l + T0T 


W h T0T 


ThThTf$T 


ZZ 


0.196 ± 0.028 


0.334 ± 0.036 


0.094 ± 0.019 


0.181 ± 0.027 


0.098 ± 0.020 


WZ 


1.058 ± 0.052 


1.087 ± 0.053 


0.447 ± 0.034 


1.006 ± 0.051 


0.248 ± 0.025 


WW 


— 


0.416 ± 0.061 


— 


0.681 ± 0.078 


0.177 ± 0.039 


it 


0.300 ± 0.057 


1.543 ± 0.128 


0.139 ± 0.038 


1.039 ± 0.105 


0.161 ± 0.041 


Z+jets 


0.112 ± 0.079 


7.3 ± 0.6 


0.168 ± 0.097 


20.3 ± 1.1 


17.9 ± 1.0 


W+jets 








37.2 ± 2.9 


6.1 ± 1.2 


case A (total) 


1.67 ± 0.11 


10.7 ± 0.7 


0.85 ± 0.11 


60.4 ± 3.1 


24.7 ± 1.6 


case B (total) 


1.49 ± 0.10 


8.0 ± 0.5 


0.76 ± 0.09 


48.4 ± 2.8 


17.9 ± 1.4 


case C (total) 


0.97 ± 0.07 


6.4 ± 0.5 


0.47 ± 0.06 


39.6 ± 2.5 


13.8 ± 1.2 


case D (total) 


0.83 ± 0.05 


4.5 ± 0.3 


0.40 ± 0.03 


31.3 ± 2.4 


9.8 ± 1.0 
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FIG. 1. Signal cross-section times branching ratio after cuts for the five channels discussed in the text: III $t (dotted), 
llThtfr (dashed), l + l + Thl^T (dot dashed), Itktk^t (dot dot dashed) and ththth I$t (solid); and for various sets of cuts: (a) 
cuts A, (b) cuts B, (c) cuts C and (d) cuts D. 
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FIG. 2. The total integrated luminosity L needed for a 3cr exclusion (solid lines) or observation of 5 signal events (dashed 
lines), as a function of the chargino mass m~+, for the three channels: III I£t (x), Uth J^t (□) and 1 + 1 + th J^t (O); and for 
various sets of cuts: (a) cuts A; (b) cuts B; (c) cuts C and (d) cuts D. 
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